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Graphene has generated a lot of research interest due to its special properties, which include 
a hydrodynamic regime. It is not yet clear however which boundary condition such a hydrody¬ 
namic current flow satisfies. The aim of this paper is to investigate the effect of different boundary 
conditions on the potential in an infinite strip of graphene, in which the electrons can be treated 
hydrodynamically. The boundary conditions on the current range continuously from no-slip to a free 
boundary. We analyse the situation for two different orientations of the source and sink, inspired by 
recent papers. We discuss which geometry is better suited for identifying the hydrodynamic regime 
and experimentally determining the boundary conditions. 


INTRODUCTION 


Correct boundary conditions are crucial but non¬ 
trivial for solving many problems in physics. One 
such problem is the flow of a fluid in a bounded 
region. For a classical fluid it is well-known that at 
the boundary between a fluid and a solid, the no-slip 
condition accurately describes the behaviour of the 
fluid [T]. The strong forces between the fluid and 
the solid make the fluid stick to the boundary, in 
all but very unusual situations. Despite being well- 
accepted nowadays, this condition was a matter of 
controversy for many years after it was introduced 
by George Stokes in the mid 19th century [2]. On the 
other hand, at the boundary between a liquid and a 
gas, the free boundary condition has to be applied 
to the liquid. 

More recently, work has been done on the hydro- 
dynamic regime in graphene, in which the electrons 
satisfy equations analogous to the fluid mechanical 
Navier-Stokes equation for an incompressible fluid 
[3] . The electrons in graphene can behave as a Fermi 
liquid and have very high mobilities. To observe this 
regime we must have a very pure sample in a spe¬ 
cific temperature range: The temperature must be 
low enough to reduce scattering of the electrons off 
the lattice, but high enough such that the mean free 
path of electrons due to collisions with each other is 
much smaller than the sample size. This can occur 
at temperatures up to room temperature [3]. The 
question remains: What conditions should be ap¬ 
plied at the boundary? Unlike in the fluidic case it 
is not clear whether to apply the no-slip condition 
here. 


Two recent papers m, m have looked for possible 
experimental signatures of the hydrodynamic regime 
in graphene. The papers both use the geometry of an 
infinite strip of width W, however different bound¬ 
ary conditions are employed. In [4] an arrangement 
is used in which source and sink are on opposite sides 
of the strip, whereas in [3] source and sink are on the 
same side. In addition [4] focusses more on the no¬ 
slip conditions at the boundary, while [3] focusses 
more on the free boundary condition. In this paper 
both geometries are investigated and the advantages 
and disadvantages of each are compared. The crite¬ 
ria are the suitability for identifying the hydrody¬ 
namic regime and the correct boundary conditions 
and for measuring the viscosity. In [4] it has been 
suggested that the negative nonlocal resistance can 
be used as an identifier of the hydrodynamic regime 
in graphene, since it does not arise for zero viscosity. 
In this article we discuss the question of how robust 
this identifier is with respect to the boundary con¬ 
ditions. 


We describe the boundary condition by a continu¬ 
ous parameter the slipping length, as in [3]. The 
possible regimes range from the no-slip condition 
/^ = 0 to the free boundary condition ^ oc. In 
this article we propose a method for experimentally 
determining the value of the slipping length. 


We focus on the potential at the boundary and 
explain its spatial behaviour and dependence on the 
relevant parameters of this problem: the slipping 
length li, and the viscosity, quantified by the 
vorticity diffusion length [3] or equivalently 
[4] . The correspondence of the dimensionless param¬ 
eters d = ^ and g = is d = Here 


V 




T] is the dynamic viscosity, e is the electron charge 
(a negative quantity), p is the ohmic resistivity and 
n is the electron number density. So g tells us the 
relative importance of viscosity and resistivity. 

In [3] and [5] it is claimed that the free bound¬ 
ary conditions are most appropriate for their par¬ 
ticular experimental set-up, since the Gurzhi effect 
is not observed. In this case the Gurzhi effect [6] 
describes the increase of the four-point conductivity 
with temperature. However, it is important to have 
a method for determining a finite slipping length 
as it is possible that the type of boundary condition 
depends on the nature of the boundary and the exact 
experimental conditions or that is large but finite, 
explaining why the Gurzhi effect is not observed. In 
the absence of conclusive theoretical arguments for 
the correct boundary conditions taking into account 
the interaction of the electrons and the boundary, 
the best that can be done is to experimentally search 
for signatures of the boundary conditions. 


V-J = 0 (1) 


+ = J (2) 

The non-linear term has been neglected, assuming 
low enough Reynolds number. The first boundary 
condition is 


Jy{x,y = ±y) = J±{x) (3) 

which defines the position of the source and sink. 
Thus we fix the current entering or exiting the strip 
though the contacts and there is no current per¬ 
pendicular to the boundary away from the contacts. 
The second boundary condition is 


SOURCE AND SINK ON SAME SIDE 


[OyJx T OxJy] 


y\y=w/2 


T 


Jx{x,y = 
Ih 


(4) 


This is the set-up used in the paper by Torre, 
Tomadin, Geim and Polini [3]. 


W 


graphene strip 

- T 


The interpretation of this last boundary condition 
is that the force exerted by the boundary of the fluid 
is proportional to the tangential velocity. is a 
phenomenological parameter of dimensions of length 
which characterizes the slipping, however it does not 
have a clear physical meaning in itself. 

In the first case we use the same set-up as in [3], 
i.e. a source and sink both at the bottom edge of 
the strip at Txq . 


We regularize via the Lorentzian with a charac- 
sink source teristic length /, so that the source and drain have 

M a finite size, as is physically sensible. Then the first 
2 Xo boundary condition becomes 


FIG. 1: geometry of the first set-up 


The origin is at the centre of the strip between 
source and sink. 


and 


J+ = 0 


(5) 


Set-up 

We use the method of [3] to solve the problem of a 
current source and sink in an infinite strip of width 
W. The governing equations are the incompressibil¬ 
ity (continuity) and Navier-Stokes equations: 


^ —e\P-\-{x — xo)‘^ P{xxo)‘^ ) 

In all the following calculations we put ^ = 1/20 
as in [4] to allow comparison with their results which 
use this value. After Fourier transforming in x, we 
can put the equations 0 to 0 in matrix form and 
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solve them as in [3] by diagonalizing the matrix. We 
find: 
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dyJx 
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y'l + ikD^y 

-i{kD^f/{l + ^-kvy/i+i/{kD,y 

1 

0 J 

_ \k\D. \ 


+<^4 


yi + (/cD,)2 

1 

0 


^kyyJl^llikD^Y 


where the Fourier transform is denoted by a hat. 

The coefficients can be determined from the 
Fourier transform of the boundary conditions ^ 
and Q. We thus obtain analytical solutions for 0, 
Jx and Jy. To obtain and Jy we Fourier trans¬ 

form numerically. 


Results 

On the following page we plot the potential at the 
lower edge = —W/2) in arbitrary units for 

p(enW) _ p{enW) _ both CaSCS WC 

plot the result for three different values of the slip¬ 
ping length, corresponding to the no-slip condition, 
a partial slip and a free boundary. When the poten¬ 
tial becomes negative, this corresponds to a negative 
nonlocal resistance. ^ 

For free boundary conditions and = \ 

(Fig. § we recover the results of the paper [3] with 
the negative resistance. Close to the source/sink the 
potential has the sign as would be expected without 
viscosity, this is due to the finite size of the elec¬ 
trode. The ideal delta-function source considered in 
[3] does not have this spike at the sources. 

The potential at infinity also has the same sign as 
in the ohmic case, the resistive term (the first term 
in (§) dominates over the viscous term (the second 
term in ([^) in the regions at large x. It is note¬ 
worthy, that the potential does not decay to zero, 
as we go to oo. Instead at positive x, for large 
viscosities, the potential tends to a positive constant 
from below, there is a flow towards the source due 
to a vortex. For small viscosities the potential ap¬ 
proaches a constant from above, as we approach the 


ohmic situation. The reason that the potential can 
tend to a finite constant as x ^ Too is that we have 
broken the left-right symmetry. So if we connected 
the ends of the strip at large x, there would be a non¬ 
zero, though small, electron flow due to the poten¬ 
tial difference since the boundary conditions define 
a sense of rotation. 

Moving away slightly from Txq the potential 
changes sign. As explained in [4] the sign change in 
the direction away from the origin is due to a vortex 
appearing. The sign-change closer to the origin can 
be thought of as due to the free boundary condition 
= 0 at the boundary y = —W/2. 


L 



FIG. 2: Illustration of the cause of negative nonlocal 
resistance 


Without any boundary and with zero potential 
the viscous term alone would create a (e.g. linearly) 
decreasing velocity profile. The boundary condi¬ 
tion = 0 changes the solution only close to the 
boundary: It lowers the current at the y = —IF/2 
boundary, as shown above. The voltage thus re¬ 
verses its sign to impose these boundary conditions. 
Keeping = \ and changing the boundary 

conditions, an additional sign change in the poten¬ 
tial appears when we decrease the slip, such that 
the sign of the voltage close to the origin is the same 
as in the Ohmic case. This is consistent with the 
previous argument. 


3 
























Geometry of [3] 



FIG. 3: 
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<t>{y = -W/2) 
= 1 or = 1 


for no-slip 
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FIG. 4: 
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(l){y = —Wl2) for partial slip ^ = 1 and 
= 1 or = 1 


(!> 



FIG. 6: 4>{y = —Wl2) for no-slip ^ = 0 and 
= 50 or ^ = 0.14 



FIG. 7: (j){y = —Wf2) for partial slip ^ = land 
= 50 or ^ = 0.14 




FIG. 5: 4>{y = —Wl2) for free boundary ^ = 1000 FIG. 8: (p{y = —lF/2) for free boundary ^ = 1000 

and = 1 or ^ = 1 and = 50 or ^ = 0.14 

p W p w 
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FIG. 9: Streamlines of the current for partial slip 
i = 1 and = 50 or ^ = 0.14 


Decreasing the viscosity such that = 50 

as in Figs. to we do not observe any negative 
resistance. We obtain something much more like the 
Ohmic case and the boundary conditions have little 
impact on the potential the reason being that we 
are exiting the hydrodynamic regime and the viscous 
term becomes negligible. 

So in this geometry we can identify two regimes 
for the dependence on the boundary conditions. For 
low viscosities (as with Dy = 50) there is hardly 
any dependence on the boundary conditions, while 
for high viscosities (as in = 1) the boundary 
conditions have an important effect. Especially the 
behaviour of the potential at the origin depends on 
the boundary conditions, with additional lobes ap¬ 
pearing for the no-slip case. This could be used to 
identify the boundary conditions. 

In the streamplot of the current (Fig. we can 
see the vortices appearing near both contact points. 
This is the cause of the negative nonlocal resistance 
which exists for these values of the parameters. How¬ 
ever the current between source an sink is much 
stronger than all other currents. 

Finally in Fig. we plot the region in the 1^- 
Djy parameter space in which (t){x^y = —IF/2) be¬ 
comes negative, i.e. the region in which we observe a 
negative nonlocal resistance. As expected, for large 
enough viscosity, around DyjW > 1 there is a nega¬ 



FIG. 10: coloured region is that region in the Ih-^v 
parameter space in which (t){x^y = —IF/2) becomes 
negative (geometry of [3]) 


tive potential regardless of the boundary conditions. 
The critical value of the viscosity for which negative 
potential is observed also decreases as we increase 
li). This agrees with our intuition that the no-slip 
condition inhibits the negative nonlocal resistance, 
since vortices at the edge of the strip would create a 
non-vanishing slip. Hence, one possible approach for 
determining the slipping length would be to find 
the critical value of D^jW at which the negative 
non-local resistance first appears. 
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FIG. 11: geometry of the second set-up 


Set-up 

In this case the governing equations are un¬ 
changed, but the boundary conditions (§-([61 must 
be replaced by 


J+ 


I I 
—e P + x‘^ 


(7) 


J- = 


I I 

—e P + x‘^ 


( 8 ) 


The problem can be solved in the same way as 
before. Due to additional symmetry, it can be solved 
by hand, as is done in the appendix. 


Results 

In the appendix we plot the potential at the 
lower edge = —IF/2) for = 30 and 

= 50. As before, in both cases we plot the 
result for three different values of the slipping length, 
corresponding to the no-slip condition, a partial slip 
and a free boundary. In the no-slip limit we expect 
the results to agree with [4]. Plotting the potential 
for different values of at /^ = 0 we find the same 


general shape of the potential. Indeed the qualita¬ 
tive picture is the same, we observe the same gen¬ 
eral form of the potential and the shape of the curve 
changes in the same way when we decrease , 

However the value of for which the curve 

r) 

adopts a particular shape is different. For example 
the critical value of at which the negative 

resistance first appears is around 47.5 in our results, 
whereas in [4] the corresponding value is 120. So 
our value of is around 2.5 times as large 

as theirs for corresponding shapes of the potential. 
This discrepancy is as yet unaccounted for. We ob¬ 
serve that if we apply the no-slip boundary condition 
to the problem, we do not necessarily find any neg¬ 
ative resistance. We also observe the potential to 
decay to zero at large x as the argument applied to 
the previous geometry does not hold here. 

This set-up allows for the experimental determi¬ 
nation of the value of the slipping length, which is a 
phenomenological parameter. A characteristic fea¬ 
ture of the graphs of the voltage in this geometry 
are the minima at Txmin, which can occur at a pos¬ 
itive or negative potential. The location of the min¬ 
imum depends on both Dj^ and so if the viscosity 
is known, this method can be used to compute the 
slipping length. Indeed, [7] find an expression for 
the viscosity of the electrons in graphene. The x- 
value of the minimum of the potential is not always 
a monotonic function of but in the case where 
it is not, a general plot of the potential along the 
boundary should allow the different possible values 
of p to be distinguished. 

We can now look at the effect that varying the 
slipping length has on the potential for the viscos- 
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FIG. 13: Streamlines of the current for partial slip 
i = 1 and = 25 or ^ = 0.2 

ity corresponding to ^ = 0.14. From Fig. l^we see 
that the minimum of the potential moves closer to 
the origin for increasing hence the vortices move 
closer to the origin. In addition, from Figs. to 
we see that the minimum value of the potential 
decreases, i.e. the negative potential becomes more 
pronounced, when we increase Z^. This makes sense 
intuitively, for no-slip conditions Z 5 = 0 we expect 
a comparatively weak vortex, since the tangential 
velocity is zero at the boundary, for free boundary 
conditions Z 5 ^ 00 this restriction does not exist, so 
the vortices become more pronounced which mani¬ 
fests itself in a stronger negative potential and the 
vortices moving closer to the origin; the latter effect 
being weaker than the other. 

Interestingly, for high enough viscosities as in 
Figs. to the dependence on the slipping 
length disappears, which is opposite to the non¬ 
hydrodynamic regime in which this happens in the 
previous geometry. Due to high viscosity strong vor¬ 
tices appear regardless of the boundary conditions. 
In fact, in contrast to the other geometry, we find 
three regimes for the dependence of the potential 
on the boundary conditions. For both high and low 
viscosities changing Z 5 has little effect. 

In Fig. [^we plot the streamlines of the current 
for parameters at which a negative nonlocal resis¬ 
tance appears. In the centre there is a strong flow 



]h_ 
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FIG. 14: coloured region is that region in the 
parameter space in which (j){x^y = —Wj2^ becomes 
negative (geometry of [4]) 

from the source to the sink. The large vortices on ei¬ 
ther side have much smaller currents flowing through 
them. Further away from the centre the flow is in the 
opposite direction, however it is significantly weaker. 

Again, in Fig. [^we plot the region in the h-Dy 
parameter space in which (t){x^y = —IF/2) becomes 
negative, i.e. the region in which we observe a neg¬ 
ative nonlocal resistance. 

As before, for large enough viscosity, around 
DytW >0.15 there is a negative potential regard¬ 
less of the boundary conditions. However the de¬ 
pendence on lb is much smaller than in the previous 
geometry. 
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Geometry of [4] 



FIG. 15: (j){y = —Wl2) for no-slip ^ = 0 and 
= 50 or ^ = 0.14 



FIG. 16: (j){y = —lF/ 2 ) for partial slip ^ = 1 and 
P(^I^ ^ 50 or ^ = 0.14 

rj W 



FIG. 18: (/)(y = —Wf2) for no-slip ^ = 0 and 
= 30 or ^ = 0.18 



FIG. 19: (/)(y = —lF/ 2 ) for partial slip ^ = land 
= 30 or ^ = 0.18 



FIG. 17: (/)(y = —lF/ 2 ) for free boundary ^ = 1000 

and = 5 Q or ^ = 0.14 

77 w 



FIG. 20 : (j){y = —lF/ 2 ) for free boundary ^ = 1000 
and = 30 or ^ = 0.18 
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CONCLUSION 


We found that the geometry of [4] leads to curves 
of a similar general shape, however it depends on the 
boundary conditions and the strength of the viscos¬ 
ity whether a negative resistance appears or not. On 
the other hand in [3] the potential can have very dif¬ 
ferent behaviour depending on the two parameters, 
with the voltage changing sign frequently along the 
boundary. When conducting experiments, the slip 
and the viscosity are a priori unknown (although 
methods exist to theoretically predict the value of 
the viscosity). The advantage of [3] is that the rich¬ 
ness of the behaviour of the voltage in the parameter 
space makes it easier to determine both parameters 
given the experimental results. If the viscosity is 
known then the geometry of [4] would make analysis 
simpler since the minimum of the potential gives an 
estimate of which regime of boundary conditions we 
are in. In fact, if we wanted to determine both the 
viscosity and the slip, we could use two markers, e.g. 
the point of the minimum together with the zero. 
An advantage of the geometry of [4] is that it gives 
a negative nonlocal resistance for smaller values of 

w • 

Another important question is whether the nega¬ 
tive nonlocal resistance is an appropriate identifier 
for the hydrodynamic regime in graphene as pro¬ 
posed in [4]. We find that the appearance of a neg¬ 
ative nonlocal resistance depends crucially on the 
size of the viscosity and the slipping length, so this 
identifier is not totally robust against changes in the 
boundary conditions. In the case of the geometry 
of [4] a more robust signature of the hydrodynamic 
regime is the presence of a minimum in the potential. 
A general trend of the influence of the slipping length 
on the negative resistance is that the negative resis¬ 
tance is smaller, or unobservable, for the no-slip con¬ 
dition compared to the others. As expected, in both 
geometries the negative resistance disappears for too 
small viscosities. However the relevant dimension¬ 
less quantity is ^ so for a small enough width W 
of the strip, the hydrodynamic regime should be ob¬ 
servable. 
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APPENDIX 


Calculation of 0 


<i){k,y) = - 


e ai{k) smh{ky) 


cro 


k‘^ 


(11) 


Due to the symmetry of the geometry in [4] this 
problem can be solved by hand. We expect the po¬ 
tential to be symmetric with respect to x and anti¬ 
symmetric with respect to y. Hence the Fourier 
transform must satisfy 


where 


ai{k) = 


I_ __ 

e 2cosh(fcVF/2) + cosh(fctV/2)f(fc) 

( 12 ) 


!i>{k,y) = ${-k,y) 


(9) 


and 


f(fc) = 


where 


4cosh(A;VF/2) -h ^ sinh(/cVF/2) 


\k\Du 




sinh(fcW/2) + cosh(fciy/2) 


(13) 


Hk,y) =-Hk,-y) ( 10 ) 


It follows that the coefficients from 0 satisfy ai = 
—a 2 and as = a 4 . We then find, using the Fourier 
transformed boundary conditions, that 


W = W 



1 

(kB,y 


(14) 


We must then Fourier transform to obtain 0, this 
is done numerically. 

In the geometry of [3] this symmetry argument 
does not work and the algebra is more tedious, hence 
we employ Mathematica to solve the problem. 
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